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I.  INTRODUCTION 


The  objective  of  this  effort  was  to  build  and  test  a  gage  capable  of 
measuring  the  specific  vertical  impulse  (peak)  imparted  by  a  high-explosive 
test  to  the  earth  or  a  berm.  The  specific  impulse  (I)  is  defined  as  the  inte¬ 
gral  of  the  pressure  (P)/time  (t)  history,  or 

I  =  fp-dt 

4 

Because  this  pressure  is  acting  on  the  soil,  it  forces  the  soil  to  move 
with  an  acceleration  (a)  equal  to 


where  mft  is  the  aerial  mass  density.  Using  this  relationship,  the  expression 
for  impulse  may  be  rewritten  as 

ft 

I(t)  =  m.  a  dt 
Jo 

Thus  for  the  case  of  constant  mass  density  at  any  time  (ti) 

I(t)  =  mA[v(ti)  -  V (to ) j 

where  t0  is  some  initial  time.  For  the  case  V(t0)  =  0,  at  any  time  during  the 
acceleration  due  to  the  high-explosive  event  (neglecting  the  material  strength 
and  assuming  a  rigid  body)  the  impulse  is  equal  to  the  momentum  of  the 
accelerated  body. 

This  is  the  principle  of  a  photopole  used  on  high-explosive  events.  A 
camera  records  the  photopole's  position/time  history.  The  velocity  of  the 
berm  is  then  inferred  from  these  data  and  the  impulse  may  be  calculated  as  a 
function  of  time.  This  technique  works  quite  well  for  a  High  Explosive  Simu¬ 
lation  Technique  (HEST)  where  the  overburden  is  directly  above  the  explosives. 
However,  in  the  Berm-Loaded  Explosive  Simulation  Technique  (BLEST),  there  are 
a  number  of  discrete  buried  charges.  The  debris  from  the  detonation  of  these 
charges  quickly  obscures  any  photopole  in  the  test  bed. 
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This  inability  to  measure  photopoles  mandated  the  development  of  some¬ 
thing  different  to  measure  impulse.  The  device  chosen  for  this  application  is 
the  time-of-flight  sphere.  It,  unlike  the  photopole,  is  used  only  to  measure 
peak  impulse.  The  peak  impulse  (Ip)  may  be  expressed,  in  the  normal  case  of 
an  initial  velocity  of  zero,  as 


I 


P 


where  Vp 


is  the  peak  or  equilibrium  launch  velocity  in  the  vertical  direction. 


The  ballistics  of  a  freely  launched  projectile  are  well-known.  In  the 
first  approximation,  neglecting  air  drag*,  simple  Newtonian  mechanics  state 
that  a  projectile  vertically  launched  in  a  1-g  downward  field  will  come  to 
rest  in  a  time  t  equal  to 


t 


u 


The  projectile  will  then  begin  to  descend  toward  the  earth  and,  if  it  lands  at 
the  same  elevation  it  took  off  from,  it  will  take  the  same  time  (td)  to 
return,  or 


2  V 

t  +  t .  =  — ^  = 
u  d  g 


where  t  is  the  total  time  of  flight.  The  above  equation  may  be  used  to 
obtain  an  expression  for  the  peak  impulse.  This  is 


The  reader  should  note,  as  will  be  shown  later,  that  air  drag  does  affect 
this  result.  However,  at  common  BLEST  launch  velocities  the  effect  is  mini¬ 
mal.  It  is  important  to  note  that  the  relationship  is  strictly  monotonic, 
i.e.,  for  any  increase  or  decrease  in  the  time  of  flight,  there  is  a  corre¬ 
sponding  increase  or  decrease  in  the  peak  impulse.  The  object  of  this  effort 
was  to  establish  that  relationship  for  earth  density  spheres  and  empirically 
verify  It.  Consequently,  a  method  for  measuring  peak  impulse  on  a  BLEST  was 
developed. 


♦Note:  The  effect  of  air  drag  is  calculated  in  Section  VI. 


II.  FABRICATION 


MECHANICAL 

The  inner  container  for  the  time-of-flight  timer  was  fabricated  from  a 
9-cm-diameter  steel  cylinder  8.5  cm  long.  One  end  had  a  screw-on  lid.  Six 
piezoelectric  crystals  were  cemented  to  the  inner  container  surface  and  con¬ 
nected  to  the  electronic  timer. 

The  outer  container  was  made  up  of  two  spun-aluminum  hemispheres  having 
a  15.24-cm  radius  and  a  1.6-mm  wall  thickness.  A  3.5-cm  aluminum  band  was 
riveted  to  one  hemisphere.  The  other  hemisphere  was  attached  to  the  band 
using  nut  plates  and  flat  socket  head  10-32NF  screws  (Figs.  1  and  2).  The  top 
hemisphere  had  a  25-mm  screw-in  plug  through  which  sand  was  poured  to  fill  the 
space  between  the  inner  container  and  the  outer  container. 
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Figure  2.  Top  hemisphere. 


III.  FABRICATION  AND  OPERATION  OF  ELECTRONICS 


TIME-OF-FLIGHT  THEORY  OF  TIMER  OPERATION— FIRST  DESIGN 

The  first  timer  built  and  tested,  unsuccessfully,  consisted  of  a  counter 
card  having  five  integrated  circuits  and  several  discrete  components.  The 
ICM7215  and  its  associated  oscillator  crystal  were  used  to  generate  all  clock¬ 
ing  signals.  While  enabled,  the  ICM7215  provided  an  oscillator  output  of 
approximately  1  kHz  to  the  4-stage  ripple-carry  binary  counter/divider  (4020). 
The  ICM7215,  as  wired,  was  configured  to  work  as  a  stopwatch.  After  being 
reset,  the  negative-going  edge  on  +  START/STOP  will  start  the  internal  counter, 
and  the  second  negative-going  edge  will  load  the  resulting  count  into  display, 
reset  the  internal  counter,  and  reinitiate  counting.  Subsequent  negative¬ 
going  transitions  would  replace  the  old  count  (i.e.,  the  current  display)  with 
a  new  "lap”  time,  reset  the  internal  counter,  and  reinitiate  counting  if  it 
were  to  occur.  The  additional  support  circuitry  surrounding  the  ICM7215  was 
designed  to  gate  only  one  start  pulse,  lock  out  further  pulses  on  the  INPUT 
line  for  a  period  determined  by  the  output  selected  from  4020,  gate  only  one 
stop  pulse,  and  lock  out  all  further  pulses  until  a  valid  RESET  occurred.  This 
timer  was  abandoned  after  bench  tests,  drop  tests,  and  explosively  propelled 
tests  showed  a  high  failure  rate  when  the  display  was  disconnected  and  recon¬ 
nected  as  is  required  during  drop  tests  and  explosively  propelled  tests. 

TIME-OF-FLIGHT  THEORY  OF  TIMER  OPERATION- FINAL  DESIGN 

The  time-of-flight  electronics  are  designed  to  meet  a  specific  need.  A 
measurement  is  made  of  a  time  increment  between  acceleration  impulses  using  an 
Armitron  Chronograph  wristwatch.  The  impulses  are  generated  by  a  positive  or 
negative  acceleration  in  any  axis  relative  to  the  device.  Piezoelectric  crys¬ 
tals  are  used  to  sense  the  presence  of  an  acceleration.  During  a  change  of 
acceleration  the  crystals  will  displace  electrons,  which  generates  a  small 
current.  The  crystals  are  arranged  in  a  three-axis  sensing  orientation  uti¬ 
lizing  two  crystals  per  axis  mounted  on  opposite  sides  of  the  canister.  Cur¬ 
rents  generated  by  the  crystals  are  fed  through  a  diode  OR  gate  to  a  voltage 
comparator.  Sensitivity  to  the  acceleration  impulses  may  be  adjusted  by 
setting  the  input  level  at  the  differential  input  of  the  voltage  comparator. 
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Output  from  the  voltage  comparator,  signaling  an  adequate  level  of  impulse,  is 
routed  to  a  flip-flop  which  is  set  at  that  time.  Output  from  the  flip-flop  is 
used  to  accomplish  two  functions.  In  the  first  function,  output  is  sent 
through  an  interface  transistor  to  the  timer  module.  Upon  receipt  of  the  sig¬ 
nal,  the  timer  initiates  the  counting  for  the  time  between  impulses.  Secondly, 
output  from  the  flip-flop  starts  the  1-s  lockout  time  to  prevent  a  false  input 
following  the  initial  pulse.  This  is  accomplished  by  an  RC  circuit  which  is 
tied  to  the  approximate  lockout  time  necessary.  The  second  acceleration 
impulse  follows  the  same  path  through  the  crystals,  the  voltage  comparator, 
and  to  the  flip-flops.  Upon  receipt  of  the  second  pulse,  the  second  flip-flop 
is  then  set.  Output  from  the  second  flip-flop  is  routed  through  an  interface 
transistor  and  then  to  the  timer.  The  second  impulse  received  sets  the  flip- 
flop  into  a  set  mode  and  will  not  allow  additional  impulses  to  pass  on  to  the 
timer.  A  complete  reset  procedure  must  be  followed  before  the  module  will 
interpret  future  impulses.  The  reset  procedure  is  a  five-step  operation  and 
is  covered  under  Time-of-Flight  Electronics  discussed  below. 

Time-of-flight  electronics— Time-of-f light  electronics  were  designed  to 
meet  qualifications  of  withstanding  a  1000-g  shock  environment  while  resolving 
the  time  of  flight  to  1/100  s.  The  electronics  are  wholly  contained  in  two 
parts.  The  first  part  contains  the  Armitron  Chronograph  wristwatch  timing  and 
interface  modules;  the  second  part  contains  the  acceleration-sensitive  switch 
and  battery  pack.  The  two  parts  are  interconnected  by  four  wire  leads.  Both 
parts  are  subsequently  installed  inside  a  heavy  metal  can  for  protection. 

The  time-of-flight  electronics  are  designed  to  be  activated  by  the  con¬ 
nection  of  the  five  interconnection  wire  leads.  Once  activated,  the  units 
should  be  able  to  operate  for  approximately  30  consecutive  days,  or  any  period 
cumulative  up  to  the  30-day  level,  which  falls  within  the  one-year  storage 
life  of  the  battery  pack.  The  development  of  the  two-part  system  allows  a  new 
battery  pack  to  be  reconnected  for  reuse  of  the  timing  module.  The  two  parts 
should  be  connected  at  a  convenient  period  before  use  and  disconnected  before 
prolonged  storage.  An  accurate  time  schedule  should  be  kept  to  determine  the 
battery  pack  usage  and  remaining  life. 
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Once  the  two  parts  have  been  connected,  the  timer  must  be  reset  by  the 
recommended  reset  procedure,  covered  later  in  this  report,  and  readied  for 
flight.  Flight  readiness  is  achieved  when  the  two  parts  have  been  joined,  the 
timer  reset,  and  both  have  been  installed  in  the  canister  along  with  the  shock 
insulating  foam  pads.  The  canister  is  then  ready  to  be  installed  in  the  sand- 
filled  sphere.  Care  should  be  exercised  to  prevent  shock  or  high  accelera¬ 
tions  being  imparted  to  the  canister  or  sphere  to  prevent  premature  starts  of 
the  timing  module.  Sensitivity  was  set  for  approximately  a  50-g  threshold  to 
start  and  stop  the  timer.  This  may  be  translated  into  a  short  drop  of  less 
than  30  cm  onto  a  hard  surface. 

The  timer  is  activated  due  to  an  acceleration  in  any  axis  which  is 
greater  than  its  threshold.  The  first  (initial)  pulse  will  start  the  timer 
running.  There  follows  a  period  of  approximately  1  s  which  is  interlocked  to 
the  timer  and  will  prevent  a  second  shock  pulse  from  stopping  the  timer.  The 
next  acceleration  impulse  following  the  interlock  time  will  stop  the  timer. 

The  timer  will  then  hold  the  time  between  impulses  on  its  display  until  exter¬ 
nally  reset.  The  timer  is  disabled  following  the  second  pulse  and  will  not 
accept  further  impulse  signals. 

To  reset  the  timer,  a  five-step  procedure,  must  be  followed.  A  30  AWG 
jumper  is  inserted  into  pin  1  and  consecutively  connected  to  pin  2,  pin  3,  pin 
4,  and  pin  5.  See  Figure  3.  At  this  point  the  display  should  read  all  zeros. 
If  not,  the  brown  wire  coming  from  below  the  watch  face  is  touched  to  the  pin 
located  at  the  lower  right-hand  corner  of  the  watch  face  and  the  procedure  is 
repeated.  The  time-of-flight  timer  is  then  ready  to  measure  the  time  incre¬ 
ment  within  the  test  environment. 

The  timer's  display  consists  of  two  lines.  The  upper  line  contains  both 
1/10-s  and  1/100-s  increments,  while  the  lower  contains  10s  and  units  of  min¬ 
utes,  followed  by  10s  and  units  of  seconds.  A  typical  time  of  flight  would  be 
contained  within  the  last  digit  of  the  lower  line  and  the  two  digits  of  the 
upper  line  of  the  display.  See  Figure  4. 


n  iiBK  m  ii4a| 


Figure  4.  Dial  of  Armitron  Chronograph  wristwatch. 
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IV.  ASSEMBLY  PROCEDURE 


The  red  wire  from  the  bottom  electronics  section  was  cut  and  taped.  Just 
prior  to  assembling  the  time-of-flight  sphere,  the  red  wire  must  be  recon¬ 
nected  by  soldering  or  using  a  twist  wire  nut  and  then  taped.  The  timer 
should  be  reset  prior  to  field  test.  To  reset  the  timer,  the  same  five-step 
procedure  discussed  previously  must  be  followed.  Pin  1  is  located  at  the 
lower  left  position  when  viewing  the  wristwatch  from  the  rear.  There  is  a 
black  spot  under  pin  1.  A  30  AWG  jumper  is  inserted  into  pin  1  and  consecu¬ 
tively  touched  to  pin  2,  pin  3,  pin  4,  and  pin  5.  See  Figure  3.  The  display 
should  now  read  all  zeros.  The  timer  is  then  ready  to  measure  the  total  time 
from  takeoff  to  landing.  The  lower  hemisphere,  which  has  the  equatorial  band 
riveted  in  place,  should  be  nearly  filled  with  sand.  The  inner  steel  con¬ 
tainer,  which  has  its  timer  reset,  is  placed  in  the  middle  of  the  sand.  Next, 

the  top  hemisphere  is  placed  inside  the  equatorial  band.  The  end  of  the  hemi¬ 

sphere  opposite  the  equatorial  band  split  is  lowered  onto  the  nut  plate,  band 
screw-hole  alignment  made,  and  then  one  10-32NF  screw  is  installed.  Because 
the  band  is  contoured  to  the  hemisphere  radius,  the  assembly  is  difficult. 

The  top  hemisphere  must  be  worked  into  the  band  by  prying  between  band  and 
hemisphere  starting  at  the  10-32NF  screw  and  working  each  side  to  the  band 
split.  See  Figure  5.  When  the  top  hemisphere  nut  plates  all  line  up  with  the 

band  screw  holes,  the  remaining  24  screws  are  installed.  The  top  plug  is  then 

removed  and  the  remaining  top  hemisphere  cavity  is  filled  with  sand  and  the 
top  filler  plug  reinstalled. 

The  assembled  spheres  are  buried  to  the  equator  for  the  field  test  shots. 
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V.  OPERATIONAL  SURVIVAL  TESTS 


The  first  timer  fabricated  consisted  of  a  counter  card  and  a  conditioner 
card  connected  to  eight  piezoelectric  crystals.  All  were  housed  in  a  steel 
inner  container  and  an  aluminum  sphere.  Drop  tests  from  the  tower  of  the  2- 
by  40-ft  shock  tube  were  conducted  at  various  elevations  from  1  to  13  m;  over 
50  percent  of  the  tests  were  failures.  The  problem  was  found  to  stem  from  the 
fact  that  the  display  unit  had  to  be  separated  from  the  timer  unit  before  the 
drop  tests  and  reconnected  after  the  drop.  This  same  timer  was  tested  for 
operation  and  survival  by  explosively  propelling  upward  the  30-cm-diameter 
sphere  containing  the  timer.  Two  explosively  propelled  tests  were  performed 
at  a  takeoff  velocity  of  15  m/s  and  25  m/s.  The  timer  failed  to  operate 
during  both  tests;  therefore,  a  new  timer  was  fabricated  using  an  Armitron 
Chronograph  wristwatch.  With  this  design,  the  watch  is  the  display  and  timer 
and  flies  with  the  sphere  as  a  unit  and  does  not  have  to  be  disconnected  and 
separated.  This  timer  system  was  dropped  from  a  height  of  13  m  over  a  hundred 
times  without  experiencing  any  failures.  Two  of  the  watch  timers  were  placed 
in  the  steel  container  but  not  the  sand-filled  sphere.  The  steel  containers 
were  then  placed  on  top  of  a  berm  which  explosively  propelled  the  timers 
upward.  Both  timers  survived  and  indicated  a  total  flight  time  of  3.85  s.  For 
the  no-drag  case,  this  gives  a  takeoff  velocity  of  19  m/s. 

After  the  above  tests,  two  watch  timers  in  two  steel  containers  were 
placed  in  one  aluminum  30-cm-diameter  sphere  with  the  void  filled  with  sand. 

The  sphere  was  placed  inside  a  38-cm-i.d.  tube  1.2  m  long.  See  Figure  6.  The 
sphere  was  then  explosively  propelled  upward.  A  photopole  and  a  1000-frame/s 
camera  were  used  to  confirm  the  takeoff  velocity  derived  from  the  time  of 
flight.  See  Figure  7.  Calculating  takeoff  velocity  from  the  time  of  flight 
for  the  no-drag  case  gave  34.7  m/s  while  the  photopole/camera  combination  indi¬ 
cated  that  the  takeoff  velocity  was  33.5  m/s.  The  error  between  the  two  time- 
of-flight  methods  was  3.33  percent. 

The  worst  error  in  determining  takeoff  velocity  using  the  camera/photopole 
is  calculated  as  follows: 

Each  photopole  alternate  black  and  white  stripe  was  50.8  mm  (2  in)  wide. 
Each  frame  of  the  1000  frame/s  camera  covered  one  stripe  due  to  the  distance 
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the  camera  was  located  from  the  photopole.  When  calculating  the  average 
velocity  as  the  sphere  passes  the  photopole,  the  maximum  probable  inaccuracy 
is  assumed  to  be  one-half  of  a  camera  frame  or  25.4  mm  as  the  sphere  enters 
the  bottom  of  the  photopole  and  +25.4  mm  leaving  the  top  of  the  photopole. 

The  photopole's  exposed  length  was  1.292  m  (48  in).  The  worst  assumed 
condition  for  accuracy  is  a  subtraction  of  50.8  mm  from  the  nominal  length  of 
the  photopole.  The  error  in  the  timing  is  presumed  to  be  plus  or  minus  the 
shutter  open  time  for  the  camera,  which  was  in  this  case  400  ys. 

Thus  the  total  decimal  error  may  be  expressed  as 

AV  _  Ax  At 
V  ‘  x 

which  is  0.037  or  3.7  percent. 
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VI.  TIME-OF-FLIGHT  AND  VELOCITY  EQUATIONS 


Up  until  now  the  subject  of  drag  has  not  been  treated.  The  travel  of 
geometric  bodies  through  air  is  obviously  affected  by  the  air.  If  they  were 
not,  sky-diving  would  have  long  ago  ceased  to  be  an  active  sport.  In  fact, 
Galileo's  original  experiment  with  free-falling  bodies,  using  the  Leaning 
Tower  of  Pisa,  is  either  only  legend  or  prevarication;  it  cannot  be  dupli¬ 
cated.  The  reason  for  this  is  air  drag. 

Baum  (Ref.  1)  presents  the  flight  equations  studied  by  Bob  Golabic  of 
RAD,  Inc,  in  Colorado  Springs.  The  following,  including  Figure  8,  are  excerpts 
from  Reference  1; 


For  the  case  going  up, 
ma  =  -D  -  W 
or 

m  ^  =  -CD  A  j  pV2  -  W  (considering  drag) 

where  m  =  mass 

CQ  *  drag  coefficient 
p  =  air  density 

A  =  cross-sectional  area  of  sphere 
W  =  weight  of  sphere 


Equation  2  may  be  rewritten  as 


1  dV  . 

g  at  “ 


1 


(1) 

(2) 


(3) 


I 


1.  Baum,  Neal,  BLEST  Impulse  Gage  Development ,  AFWL-TR-80-13,  Air  Force 
Weapons  Laboratory,  Kirtland  Air  Force  Base,  New  Mexico,  October  1980. 


r 


D  =  drag 
W  =  weight 

ma  =  mass  x  acceleration 

Figure  8.  Forces  acting  on  vertically  launched  ball. 
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or 


V 

f 


dV 


pCPA 

"2W 


V2  +  1 


t 


;o 


gdt 


where  V0  is  the  launch  velocity.  If  one  defines  Co  as 

pCqA  3CqP 
Co  =  “2Wg  =  Sp^gr 

and 


Ci  =  arc  tan  (VoCo) 


(4) 

j 


(5) 

(6) 


then,  assuming  a  constant  drag  coefficient, 
-gt  *  jU-  £arc  tan  (VoC0  -  Ci)J 


(7) 


To  determine  the  total  ascent  time  (Tu),  it  is  noted  that  V-K). 


Thus, 


lu  ■  §tg  ‘  gt  arc  U"  <V'C,) 


(8) 


To  determine  the  descent  time  (t^),  the  peak  height  (Z)  of  the  ascent 
must  be  determined.  Thus,  rearranging  Equation  7,  it  can  be  shown  that 

3F  a  h  tan  ^Cl  ‘  C°gt) 


(9) 


or,  when  integrating  the  right  side  between  0  and  t 

Z  =  -  In  (cos  Ci)  (10) 

Now,  in  the  case  of  a  downward  trajectory 

ma  =  -D  +  W  (11 ) 
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or 


dV 


1  -  C0"V 


7yT  =  gdt 


if 


then 


k  =  r* 

Lo 


k2dV 

k2  -  V'2  “ 


gdt 


or 


v  = k 


e¥- 1 


1  +  e 


This  equation  may  be  rewritten  in  the  form 

dt 


L. 

* 

dx  =  k 

f  d 

tan  h 

e 

f 

^  / 

0 

0 

or 


Co2g 


In  cos  h  Cogtd| 


Returning  to  Equation  10  and  substituting  the  value  of  Z, 
td  =  arc  cos  h  j^cos  [arc  tan  (V0C0)])  | 


2V 


*u  +  V 


The  equation  for  the  no-drag  case  is  t  *  =^-  where 
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(12) 

(13) 

(14) 

(15) 

(16) 

07) 

one  finds  that 

(18) 

t  =  total  flight  time 


1 


Appendix  A  shows  a  sample  calculation  computing  time  of  flight  and  velo¬ 
city  (Ref.  2).  Table  A-l  in  the  appendix  lists  time  of  flight  versus  launch 
velocity  for  an  earth  density  sphere.  Figure  9  is  a  graphical  plot  of  time 
of  flight  versus  launch  velocity  with  atmospheric  drag. 

2.  Horner,  I.  B.,  Fluid  Dynamic  Drag ,  published  by  the  author,  1958. 


VII.  CONCLUSION  AND  RECOMMENDATIONS 

The  first  conclusion  has  to  do  with  accuracy.  The  experiments  in 
launching  the  sphere  indicated  that  the  velocity  inferred  from  the  time  of 
flight  was  accurate  to  3.6  percent.  This  is  somewhat  difficult  to  interpret 
because  the  accuracy  of  the  velocity  measurement  (photographic)  was  only 

3.7  percent.  The  time  of  flight  was  at  least  as  accurate  as  the  standard  used 
to  measure  it. 

Other  measurements  were  made  using  the  timer  only  and  not  the  sphere*. 
These  were  done  on  a  shaped  HEST,  part  of  a  test  called  HYBD-1  (21).  These 
measurements  are  summarized  in  Table  1  from  Reference  3.  Here  the  average 
agreement  with  the  time-of-fl ight  and  other  measurements  of  impulse  is 

2.8  percent.  This  is  better  than  the  normal  measurements  of  pressure,  accel¬ 
eration,  and  photopoles. 

TABLE  1.  TOTAL  IMPULSE  AS  A  FUNCTION  OF  MEASUREMENT  TYPE 
AND  ZONE  OF  THE  SHAPED  HEST 


Total  Impulse,  MPa-s 


Measurement  type 

Time  of  Flight 

Pressure 

Acceleration  1 

Photopole 

Zone  8 

0.0135 

0.0136 

0.0134 

0.0136 

Zone  5 

0.0126 

0.01245 

0.0129 

— 

Zone  6 

0.0108 

0.0119 

— 

— 

(Ref.  3) 


The  theoretical  accuracy  of  the  device  can  be  simply  determined  from  the 
accuracy  of  the  timer.  This  is  ±0.01  s.  For  an  average  time  of  flight  (7  s), 
the  accuracy  in  velocity  is  0.14  percent.  This  presumes  that  all  the  constants 
usej  in  correcting  the  drag  are  correct.  This  is  obviously  not  true. 


3.  Palmer,  D.  G.,  HYBD-1  (21)  Final  Data  Report,  NMERI  TA9-14,  Air  Force 
Weapons  Laboratory,  Kirtland  Air  Force  Base,  NM,  December  1981. 

*While  it  is  noted  that  air  drag  is  a  minimal  effect  at  these  velocities,  the 
author  still  recommends  the  use  of  a  sphere.  This  is  for  the  simple  reason 
that  the  drag  effects  on  a  sphere  are  known. 
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The  conclusion  is  that  the  instrument  being  tested  is  more  accurate  than 
the  standards  used  to  test  it.  Since  the  worst  velocity  inaccuracy  observed 
was  of  the  order  of  4  percent,  it  can  be  stated  that  the  time-of-f light  meas¬ 
urement  of  the  velocity  is  at  least  this  accurate.  Thus  Figure  9  may  be  used 
by  the  experimenter  to  determine  launch  velocity. 

Table  1  also  indicates  another  advantage  of  the  time-of-fl ight  gage.  The 
reason  some  of  the  measurements  were  not  given  is  because  of  failure.  In  the 
case  of  the  photopoles,  failure  resulted  from  obscuration.  In  the  case  of 
acceleration,  cable  failure  was  the  cause.  There  were  no  failures  in  time-of- 
flight  measurements.  Thus  Table  1  indicates  relative  reliability  as  well  as 
accuracy. 

The  accuracy,  reliability,  and  ease  of  installation  indicate  to  the 
author  that  the  time-of-fl ight  gage  should  be  used  in  all  tests  where  appro¬ 
priate.  This  includes  both  BLEST  and  HEST. 


24 


APPENDIX  A 

TIME-0F-FL1GHT  AND  VELOCITY  SAMPLE  CALCULATIONS 

Determine  drag  coefficient,  Cd,  from  1  m/s  to  160  m/s.  Calculate  Rey¬ 
nolds  number,  R,  from  1  m/s  to  160  m/s  using  information  from  Reference  2. 

R  =  ^i_ 

r 


where 

V  =  ft/s,  velocity 
1  =  body  length,  1  ft 

/~  =  kinetic  viscosity  =  1.56  x  10“ 4  ft2/s 


Table  A-l  from  Reference  2  shows  a  constant  CQ  = 
and  4  x  105.  At  4  x  10s,  C^  drops  to  0.1  and  remains 


0.47  between  R  =  1 0 3 
constant  beyond  10®. 


Determine  velocity  at  which  R  =  4  x  105  is  the  crossover  point  from 
CD  =  0.47  to  Cn  *  0.1. 


R  =  — 

r 


r  r 

i 


V  =  (400,000)  (0.000156)  =  62.4  ft/s 

or  =  19.02  m/s 

R^R°sea  level  =  anc*  vei°cit.y  is  directly  proportional  to  R.  The  crossover 
velocity  is  V  *  (0.9) (19. 02)  =  17  m/s  at  5000- ft  altitude.  Therefore, 

CQ  =  0.47  should  be  used  for  all  velocities  less  than  17  m/s  and  C^  =  0.1  for 
all  velocities  greater  than  17  m/s. 

The  following  is  a  sample  calculation  determining  time  up,  t  ,  time  down, 
td,  and  total  time  of  flight  for  velocity  V0  =  100  m/s  taking  drag  into 
account. 
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TABLE  A-l .  TIME-OF-FLIGHT  VERSUS  LAUNCH 

VELOCITY  FOR  EARTH  DENSITY  SPHERES 


Drag 

coefficient, 

CD 

Sphere 
liftoff 
velocity, 
V0  m/s 

\ip 

tdown 

Total 

Wag 

Total 

tno-drag 

%  difference 

‘drag  versus 
^no-drag 

0.47 

4 

0.41 

0.56 

0.97 

0.82 

15.46 

0.47 

8 

0.82 

0.79 

1.61 

1.63 

1.24 

0.47 

12 

1.22 

1.25 

2.47 

2.45 

0.81 

0.47 

16 

1.63 

1.67 

3.30 

3.27 

0.91 

0.1 

20 

2.04 

2.04 

4.08 

4.08 

0 

0.1 

40 

4.08 

4.08 

8.16 

8.16 

0 

0.1 

60 

6.10 

6.11 

12.21 

12.24 

0.25 

0.1 

80 

8.11 

8.14 

16.25 

16.33 

0.49 

0.1 

100 

10.10 

10.16 

20.26 

20.41 

0.73 

0.1 

120 

12.00 

12.08 

24.08 

24.49 

1.70 

0.1 

140 

13.96 

14.10 

28.06 

28.57 

1.82 

0.1 

160 

15.86 

16.06 

31.92 

32.65 

2.29 

Note:  Values  in  seconds  except  where  indicated. 
Air  density,  p  =  1.22  kg/m3 
Soil  density,  ps  =  1760  ks/m3 

Gravitational  acceleration,  g  =  9.8  m/s2 
t  =  tu  +  tj  (time  up  and  time  down) 

Sphere  radius,  r  =  0.1524  m  (6  in) 

t  = 

no-drag  g 


r 
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From  Tabl e  fl-1,  where 

V  ,  -  100  m/s 


air  dens i  t.y  .=1.22  kg/m  ' 

soil  density  p  =  1760  kg/m 

gravitational  acceleration  g  =  9.8  m/s' 
sphere  radius  r  =  0.1524  m  (6  in) 


„  3tD  ■  (3) (0. 1 ) ( 1 . 22) 

L,J  '  8~gr  '  (8)  ( 1  760)  ( 9.8)  (0. 1  524  ) 


0.0017405 


Time  up,  t  =  l/gC}  arc  tan  (V,C0)  from  Equation  8. 

lu  =  1978) (0.001 7405)  arc  tan  [H00)(0. 001  1  7405)] 
t  -  ( 58. 6273) (0. 1 723237 )  =  10.102876  s 
Time  down,  td  -■  1/gC-  arc  cosh  cos  [arc  tan(V  C  )J  from  Equation  18. 

ld  =  (9. 3) CO. 0017405)  arC  C0Sh  O-0150337) 

t.  =  58.673  arc  cosh  1.0150337 
d 

cosh'1  1.0150337  =  0.173183 


Td  --  (53. 673)  (0. 1 731 83 )  =  10.161168  s 


Total  flight  time,  t  =  t  +  t  ,  =  10.10  +  10.16  =  20,26  s 
J  u  d 

For  the  no-drag  case,  the  total  flight  time  is  t  =  2V  /  g  ••  ( 2 )  ( 1 00 ) /9 . 8  = 
20.41  s. 

If  the  no-drag  total  flight  time  equalled  the  drag  case  total  flight 
time,  the  takeoff  velocity  V :  would  be 


V-  =  =  (2Q-2|H?--8-i  =  99.274  m/s 

The  percent  error  using  the  no-drag  case  versus  the  drag  case  would  be 
100  -  99.274  x  100 


error 


100 


=  0.726  percent 
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